Magnetic phase diagram of the antiferromagnetic pyrochlore Gd2Ti207 
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Gd2Ti207 is a highly frustrated antiferromagnet on a pyrochlore lattice, where apart from the 
Heisenberg exchange the spins also interact via dipole-dipole forces. We report on low-temperature 
specific heat measurements performed on single crystals of Gd2Ti207 for three different directions of 
an applied magnetic field. The measurements reveal the strongly anisotropic behaviour of Gd2Ti207 
in a magnetic field despite the apparent absence of a significant single-ion anisotropy for Gd^^ . The 
H — T phase diagrams are constructed for ii" || [111], H || [110] and H \\ [112]. The results indicate 
that further theoretical work beyond a simple mean-field model is required. 

PACS numbers: 75.50.Ee, 75.30.Gw, 75.30.Kz, 75.60.Ej 



Amongst the members of the titanium pyrochlore ox- 
ide family with the general formula R2Ti207, the mag- 
netic properties of which are all presently the subject of 
intense investigations, gadolinium titanium oxide (GTO) 
holds a unique position. In contrast to the other members 
of this series of compounds, where magnetic anisotropy 
plays a crucial role and results in the appearance of some 
most unusual ground state configurations, including the 
"spin-ice" state GTO is the only compound where 
the magnetic rare earth ions have a negligible single-ion 
anisotropy. This is because in Gd2Ti207, the Gd'^^ ions 
are in a state with S — 7 /2 and L — Q. The magnetic 
ions are antiferromagnetically coupled. Therefore, apart 
from the ever present dipole-dipole interactions, GTO 
may be regarded as a realisation of an ideal Heisenberg 
antiferromagnet on a frustrated pyrochlore lattice. 

The considerable theoretical interest in the Sd py- 
rochlore model and its 2d analogue (a checkerboard lat- 
tice), is mostly due to the fact that these models, both in 
classical and quantum descriptions, do not develop con- 
ventional magnetic order down to T = K 2]. Their 
ground state is instead described as a degenerate spin- 
liquid Whether the weaker perturbations such as a 
single ion anisotropy, dipolar interactions 0, 0], small 
lattice distortions^a , further- neighbor interactions |^ or 
a magnetic field |8j break the degeneracy and force the 
system to adopt a particular ground state configuration 
is the subject of current investigations. In this context 
it is extremely important to identify different magnetic 
phases in real magnetic materials and to establish the 
phase diagrams. 

Another unusual aspect of the thermodynamics of frus- 
trated magnets, which may have direct consequences for 
GTO, has been pointed out recently by Zhitomirsky 0. 
Because a significant enhancement of the magnetocaloric 
effect is expected for frustrated magnets near the sat- 
uration field, GTO may have potential as an adiabatic 
demagnetisation refrigerator. An enhancement of the 
magnetocaloric effect in model systems is attributed to 
the presence of a significant number of so-called "soft 
modes" . The existance of such soft (and zero) modes 



in a magnetic material such as Gd2Ti207, still needs to 
be established. To date, the only direct experimental 
observation of soft modes has been made with inelas- 
tic neutron scattering in another frustrated compound 
ZnCr204 [13 ■ 

In this letter, we present the results of specific heat 
measurements on Gd2Ti207 in an applied magnetic field. 
By using single crystal samples we have been able to 
prove that there is an unusually high degree of anisotropy 
built into the ground state and to construct the H — T 
phase diagrams for three different directions of magnetic 
field, H\\ [111], H\\ [110] and II [112]. 

The single crystals of Gd2Ti207 were grown by the 
floating zone technique, using an infrared image fur- 
nace [ll|. The principal axes of the samples were de- 
termined using X-ray diffraction Laue photographs; the 
crystals were aligned to within an accuracy of 2-4°. The 
specific heat measurements were performed in the tem- 
perature range 0.39 to 5.0 K in a field of up to 9 T using 
a Quantum Design PPMS calorimeter equipped with a 
■^He option. In order to estimate the lattice contribution 
to the specific heat we have measured the specific heat 
of a single crystal of Y2Ti207, a nonmagnetic compound 
isostructural to GTO. For all the low- temperature mea- 
surements reported in this letter (T < 5 K) the lattice 
contribution was found to be insignificant. 

We have measured the specific heat both as a function 
of temperature in a constant magnetic field and as a func- 
tion of applied field at constant temperature. Because of 
the relatively high values of the specific heat of GTO and 
its poor thermal conductivity at low temperatures, the 
measurements were restricted to small (< 1 mg) plate- 
like samples with a typical thickness of 0.1 mm. For some 
orientations of applied field, the demagnetizing factors 
were significant and had to be taken into account |l2l |. 
This has been done by measuring the magnetization of 
the same or similarly shaped samples with either a vi- 
brating sample or a SQUID magnetometer. 

The temperature dependence of the heat capacity di- 
vided by temperature, C{T)/T, measured in zero field is 
shown in Fig.^ Two sharp peaks of nearly the same am- 
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FIG. 1: Temperature dependence of the specific heat divided 
by temperature, C(T)/T, measured on the single crystal of 
the Gd2Ti2 07 in a zero applied field. The inset shows the 
temperature dependence of the magnetic entropy. 
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FIG. 2: Temperature dependence of the specific heat, C(T), 
of Gd2Ti207 at various values of an external magnetic field 
applied along the three different crystallographic directions, 
_H"||[flf], //||[112] and //|t[110]. 



plitude observed at TAri=1.02(2) K and TAr2=0.74(2) K 
mark the temperature of two consecutive phase transi- 
tions. In Fig. 1 two data sets are shown in order to 
demonstrate the high degree of reproducibility in these 
measurements. These data should be compared with the 
results of previous measurements performed on powder 
samples Raju et al. 0| reported a broad peak 

centred around 2 K (absent in our data) and a sharp 
peak just below 1 K, which were attributed to the pres- 
ence of short range correlations and the development of 
long range magnetic order respectively. More recently, 
however, Ramirez et al. [l^ have reported the presence 
of another sharp peak at around 0.7 K. While our data 
are similar to those presented in Ref . ,14j , they differ sig- 
nificantly in the absolute values of C(T) from the values 
reported in Ref. |l3|. For example, at the temperature 
of the upper phase transition, rjvi = 1.02 K in our data, 
the specific heat reaches a value of 22 Jmol^^K^^, while 
the values reported in Ref. j.l3i] are approximately half 
this quantity. Moreover, Bonville et al. [isj have very 
recently repeated the measurements of specific heat on a 
powder sample of GTO and reported values which are in 
a good agreement with our data. 

In terms of the overall magnetic entropy, Raju et al. 
observed a significant magnetic specific heat up to about 
30 K, while Ramirez et al. have recovered 92% of the 
expected entropy, 2i?ln(25' + 1) = 34.6 Jmor^K^i, be- 
tween and 5 K. The differences in the specific heat 
data, as well as in the ac magnetic susceptibility be- 
haviour, have been suggested to arise from the extreme 
sensitivity of the low-temperature properties to the sam- 
ple preparation technique [l^ . The inset in Fig. 1 shows 
the temperature dependence of the magnetic entropy, ob- 
tained by extrapolating the C{T)/T to K and numer- 
ically integrating it versus temperature. Approximately 



90% of the expected value is recovered when the tem- 
perature reaches 5 K. Additional high temperature mea- 
surements have shown that the entropy is almost fully 
recovered by 10 K, a temperature which corresponds to 
the reported Curie- Weiss temperature of GTO At 
that point the magnetic contribution to the specific heat 
becomes negligible when compared to the lattice contri- 
bution. The fact that a majority of entropy is recovered 
only at T >Tni signifies the importance of short-range 
correlations in GTO. 

Powder neutron diffraction measurements performed 
on GTO sample by Champion et al. in zero field |l6l |. 
have revealed a complex noncoUinear structure with a 
magnetic propagation vector ofk = (51^)- According 
to these measurements, the GTO magnetic structure (the 
so-called P state) consists of kagome planes with "g = 0"- 
type order with interstitial sites carrying zero magnetic 
moment; GTO has been suggested to be only partially 
ordered even at 50 mK The observed recovery of 

the full magnetic entropy conflicts with this suggestion. 
On the other hand, mean-field theory Q predicts that an 
ordered four sublattice structure (F state) competes with 
the P state and becomes more stable as the temperature 
is reduced. Therefore we suggest that at Ti^i GTO un- 
dergoes a transition from a paramagnetic to the P state, 
while at Tn2 the transition is to the F state. 

Our lowest experimentally accessible temperature of 
0.39 K is about 50% of Tn2- We cannot therefore re- 
liably extrapolate the C{T) dependence down to K. 
From the data shown in Fig. 1, however, it is quite clear 
that between 0.39 K and Tn2, C{T)/T grows approx- 
imately linearly with temperature, which implies a 
dependence for the low-temperature part of C{T). A 
similar dependence has been observed in another mag- 
netically frustrated system, SrCr9pGai2-9pOi9 (SCGO) 
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FIG. 3: Field dependence of the Gd2Ti207 specific lieat di- 
vided by temperature, C{T)/T, for tliree different directions 
of an applied magnetic field. 



with a kagome lattice [l^. 

The temperature and field dependencies of the spe- 
cific heat, measured on the single crystals of GTO for 
the three different directions of applied magnetic field, 
H II [111], H II [112] and H \\ [110] are shown in Figures 
2 and 3. The application of a small external field causes 
the convergence of the transition temperatures, T/vi and 
Tn2, regardless of the field orientation. In a field of about 
1.5 T, the two transitions merge as seen on the left- 
hand panels of Fig. 2. In a small field and low-T, C(T) 
still follows a dependence; above 1.5 T, however, there 
is no dependence in C{T). This is in contrast with 
SCGO, whose specific heat shows almost no field depen- 
dence 113. For fields above 1.5 T, significant differences 
in the positions of the peaks in the C{T) curves for GTO 
are clearly visible (see right-hand panels of Fig. 2). 

The field dependence of the specific heat shown in 
Fig. 3 has been measured in order to determine pre- 
cisely the positions of the phase transitions lines at low 
T, where in high fields the specific heat is nearly temper- 
ature independent. The higher-field peaks observed in 
these C{H) curves at around 5 to 6 T, mark a transition 
to a saturated (paramagnetic) phase. The low-T peaks 
found at approximately half the saturation field corre- 
spond to an additional field-induced transition discussed 
below. The most intense peaks in C{H) are found in 
the higher-T low-field region. They reach values above 
70 Jmol~^K~^ and correspond to a transition into an in- 
termediate temperature phase. It is interesting to com- 
pare this C{H) data with the results of Monte Carlo sim- 
ulations for classical spins on the pyrochlore lattice near 
the saturation field where the presence of a macro- 
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FIG. 4; Magnetic phase diagrams of the Gd2Ti207 for three 
different orientations of an applied magnetic field. 



scopic number of soft modes has been predicted. This 
comparison is not straightforward, as in the simulations 
only the exchange interactions have been included, which 
precludes the system from entering a state with long 
range magnetic order Q . Fig. 3 shows that for H > Hgat 
the measured field dependence of the specific heat does 
not follow a simple exponential decay expected for a non- 
frustrated antiferromagnct. Instead it remains almost 
field-independent for H < 7 T and then falls rapidly at 
higher fields. This observation may suggest a large in- 
fluence of the soft modes on the specific heat of GTO 
around the saturation field. 

No signs of history dependence have been detected 
throughout the course of our measurements, e.g. the 
values of the measured specific heat were identical for 
increasing/decreasing field and temperature. This sug- 
gests a second order nature for all the phase transitions 
in GTO in contrast with the observations of Bonvillc et 
at 15], which imply the upper phase transition at T^i 
has a weak first order character . 

The magnetic i? — T phase diagrams for GTO obtained 
by plotting the temperature of various phase transitions 
observed in both the C vs. H and C vs. T measurements 
are shown in Fig. 4. By analysing these phase diagrams, 
the highly anisotropic nature of the magnetic behaviour 
of GTO becomes more apparent. To start with, when 
H II [111], the magnetic field favors the onset of mag- 
netic order in GTO. For example, in a field of about 
4 T the transition temperature is 1.14(2) K compared to 
1.02(2) K in zero field, while the transition to a saturated 
phase measured at 0.57 K occurs at about 5.7 T. When 
H\\ [112], an external field inhibits the onset of long range 
magnetic order: at 4 T the transition is reduced to 0.82 K 
and the saturation field at 0.57 K does not exceed 5.2 T. 
Application of the field along the [110] direction has an 
intermediate effect. 

Even more remarkably, an additional phase transition 
induced by the external field at about half of the sat- 
uration field is clearly present at low temperatures for 
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H II [111], but is absent for H \\ [112]. For H \\ [110], this 
phase transition is present, but it becomes less and less 
pronounced as the temperature decreases (see bottom 
panel of Fig. 3). The origin of this phase transition is not 
clear at present, although the fact that it happens at a 
half of Hsat strongly suggests that it is likely to be related 
to a collinear spin state, where three spins on each tetra- 
hedra are aligned with the field and the fourth is pointing 
in the opposite direction. A stabilisation of the collinear 
phase by quantum and thermal fluctuations has been pre- 
dicted to take place at H = Hsat/2 for Heisenberg an- 
tiferromagnets on pyrochlore and frustrated square lat- 
tices, as well as at H = Hsat/'^ for kagome and garnet 
lattices 0. Whether or not GTO, whose magnetic in- 
teractions are more complex than a simple Heisenberg 
model, does indeed adopt a collinear order in a magnetic 
field remains to be established with more direct probes of 
the magnetic structure, such as neutron diffraction mea- 
surements. 

When compared to the phase diagram obtained from 
the C{H,T) measurements on a powder sample [l^ . 
Fig. 4 demonstrates the advantages of using single crys- 
tals. While it has been claimed that there are as many 
as four different ordered phases in GTO 0], our mea- 
surements show that for any given direction of magnetic 
field there are at most three ordered phases. The phase 
diagram reported by Ramirez et al. |l4| seems to consist 
of a superposition of the single crystal phase diagrams. 

When compared to the theoretical predictions based 
on a simple mean-field model 0|, our phase diagrams 
reveal significant discrepancies. For example, for the field 
along [111], theory predicts only one transition, while 
our measurements clearly show two. These observations 
strongly suggest the need for further development of the 
theoretical model, which most likely should include some 
degree of magnetic anisotropy. 

Recently Hassan et al. [13 have reported the presence 
of a surprisingly large anisotropy in GTO below a tem- 
perature of 80 K, which they observed by measuring the 
ESR spectra. The anisotropy reveals itself as a well pro- 
nounced shift in the position of one of the resonance lines 
from its maximum value for H \\ [111] to its minimum for 
HI. [111]. The shift amoimts to 3.6 T when measured 
at a frequency of 54 GHz. No change in the position of 
the resonance lines was observed when the field was ro- 
tated in a plane perpendicular to the [111] direction p^ . 
The same authors have also reported dc-torque measure- 
ments [2oj . which have shown a gradual development of 
macroscopic anisotropy below 80 K. 

Motivated by these results, we have measured the tem- 
perature dependence of the magnetic susceptibility of a 
GTO crystal for H \\ [111] and _L [111] using a SQUID 
magnetometer. No significant anisotropy in the suscep- 
tibility has been detected At present, neither the 
origin of the anisotropy observed in the ESR measure- 
ments [19I ]. nor its presence at a relatively high temper- 



ature have a satisfactory explanation. Any explanation 
should obviously take into consideration more than a sim- 
ple Heisenberg exchange interaction between the nearest 
neighbor spins. 
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